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Abstract 

This work considers a Blasius boundary value problem with inhomogeneous 
lower boundary conditions f(0) = 0 and f'(0) = - X with X strictly 
positive. The Crocco variable formulation of this problem has a key term 
which changes sign in the interval of interest. It is shown that solutions of 
the boundary value problem do not exist for values of X larger than a 
positive critical value X . The existence of solutions is proved for 
0 < X < X by considering an equivalent initial value problem. However, for 
0 < X < X , solutions of the boundary value problem are found to be 

nonunique. Physically, this non-uniqueness is related to multiple values of 
the skin friction. 
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Introduction 


The simplest example of Prandtl's boundary layer theory is the flow along 
a thin, semi-infinite flat plate (y = 0, x > 0) immersed in an 
incompressible liquid of low viscosity flowing with a given constant velocity 
in the x direction (Schlichting, 1955)* The boundary layer flow in this 
situation is governed by the Blasius equation 

+ ff~ = 0 (1) 


subject to the boundary conditions 

f (0) = 0, f'(0) = 0, and f'(«) = 1 (2) 

where the dimensionless streamfunction f(n) is ^(x,y)/(2vx) ^ with 
being the dimensional streamfunction, v the kinematic viscosity, and B the 
dimensionless similarity variable y/(2vx) ^ . The classical solution of 
Blasius of (1) and (2) was based on a coordinate perturbation method. It 
involved a power series expansion in n about the origin and an asymptotic 
expansion for large n. For the coordinate perturbation method to be 
consistent, the regions in which the inner and outer expansions are valid must 
overlap so the expansions can be matched. A rigorous analysis of the Blasius 
problem (as well as the more general Falkner-Skan problem) was not available 
until Weyl (1941) used function theoretic techniques to study the existence 
and uniqueness of solutions of the Blasius equation. Weyl showed that the 
radius of convergence of the series solution of the Blasius equation lay 
between 2.08 and 3.11. This approach has been further extended by Hussaini 
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(1971) to treat the case of the Blasius equation with three-point boundary 
conditions, a type of situation which arises in connection with the mixing 
problem behind an infinite step. 

A simple proof of Weyl's results for the Falkner-Skan equation has been 
provided by Coppel (1960) using the theory of differential equations. Coppel 
approached the Falkner-Skan equation 

f'" + ff~ + 3(1 - f'h = o O) 

by embedding the boundary conditions at y = 0 in (2) in the class of 

problems with boundary conditions f(0) = a and f'(0) = b where a and b 
are arbitrary non-negative constants. He obtained bounds on f — (0) of the 
form 

y e < [f~(o )] 2 < y 3 + 1 . ( 4 ) 

The special cases of the Blasius equation (3=0) and the Homann equation 
(3 = 1/2) were considered separately. 

For the Blasius problem, the value of f'"(0) is .469600... . It should 
be noted that once this value has been determined, the Blasius two-point 
boundary value problem can be solved in a straightforward manner as an initial 
value problem. For the Blasius problem, f"(0) can easily be obtained using 
homotopy methods (Rosenhead, 1963). The success of the homotopy technique 
depends crtically, however, on the identically zero values of f(0) and 
f'(0). 

Callegari and Friedman (1968) have noted that while Weyl's iterative 
technique for solution of the Blasius problem is convergent, in practice it 
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cannot be carried beyond two iterations. As an alternative, they considered 
the problem in Crocco variables. In this approach, the tangential velocity 
u = f'(ri) becomes the new independent variable while the new dependent 
variable is the shear stress g(u) = The Blasius problem now belongs 
to a class of nonlinear boundary value problems of the form 



g(u)g"(u) + 2uk(u) = 0, 

0 < u < 1 

(5) 


g'<0) = 0 


(6a) 


g(l) - o 


(6b) 

where k(u) is 

a non-negative continuous 

function for u 

in the closed 

interval [0,1]. 

For the Blasius problem, 

k(u) = 1. Using 

techniques of 


analytic function theory, Callegari and Friedman prove the existence, 
uniqueness, and analyticity of solutions of the Blasius problem. Their proofs 
are critically dependent, however, on the fact that the function uk(u) in 
their singular differential equation does not change sign on the interval 0 < 
u < 1. Thus, if k(u) is non-negative and the solution exists at a point ug, 
then g'(ug) <0 so g(u) is a decreasing function. 

The analyticity of solutions of equation (5) with k(u) = 1 and the two 
point boundary conditions 

g(«!) - 0, g'<u 2 ) = 0 (7) 

has also been discussed by Callegari and Nachman (1978). This problem has 
applications to the flow behind weak expansion and shock waves and the 
boundary layer flow on a conveyor belt. 
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The present work considers a two-point boundary value problem consisting 
of the Blasius equation (1) and the nonhomogeneous boundary conditions 

f(0) = 0, (8) 

f'(0) = - X, 

together with the outer boundary condition 

f '(“) = 1 (9) 

where X > 0. The Crocco variable representation of this boundary value 
problem is 

gg~ + (x - X) = 0, 0 < x < 1 + X (10) 

g'(0) = 0 and g(l + X) = 0 (11) 

where x = {f'(ri) + X} and g(x) = f"(n). Unlike the usual Blasius problem 
with homogeneous conditions at n = 0 which leads to (5), the second terra in 
equation (10) changes sign in the interval of interest. Hence, g(x) will not 
be monotone on the interval. A critical condition present in previous work on 
the Blasius problem is thus violated here. 

In the next section, we show that if a solution to the boundary value 

problem (10) and (11) exists, then there is an upper bound X such that 

X < X • Section 3 proves the existence of a continuous solution of (10) and 

(11) when 0 < X < X • However, as shown in section 4, a solution for fixed X 
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in this range is not unique. The existence of X and the non-uniqueness of 
solutions of the boundary value problem are previously unsuspected properties 
of (10) and (11) directly traceable to the fact that x - X is not purely 
non-negative on the entire interval of interest. 


2. Dependence on the Parameter X 

Suppose that a solution g(x) exists for the boundary value problem (10) 
and (11). Then, the identity 

igg')' - gg" + ( g ') 2 

can be used to write equation (10) in the form 

(gg')' + (X - X) = (g')2. (12) 

Let g(0) = a. Then, using the fact that g'(0) = 0, integrating equation 
(12) from 0 to x and integrating the result a second time from 0 to 
X = 1 + X gives that 

X 2 (X - 3X)/6 = <x 2 /2 + / dx / g' 2 (s)ds. (13) 

0 0 

The right-hand side of equation (13) is intrinsically positive and thus 


X 2 (X - 3X) > 0. 
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Setting X = 1 + X now gives that 


X < V2 • (14) 

Thus, if a solution to the boundary value problem (10) and (11) exists, X 
must be bounded above by a finite number X , and an upper bound on X is 
given by (14). In fact, while (14) shows that a bound on X is required, it 
is a rather crude upper bound as it is obtained by simply noting that the 
right-hand side of (13) is positive. A sharper upper bound can be obtained by 
an analysis which obtains a sharper positive lower bound for the right-hand 
side of (13). This analysis will be done in a subsequent paper. In section 4 
of the present paper, numerical results are given which show that 
X* = 0.3541078... . 


3. Existence of Solutions 

Consider now the existence of solutions of the boundary value problem (10) 
and (11) for 0 < X < X . The required result will be shown by proving the 
following theorem for an equivalent initial value problem. 

Theorem: There exists at least one positive constant a such that the 

initial value problem consisting of equation (10) and the initial conditions 

g(0) = a, g'(0) = 0 (15) 


has a continuous solution which is positive for 0 < x < 1 + X and 
g(l + X) = 0. 
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Proof: As a preliminary step, the initial value on g in (15) will be 
normalized to unity and dependence on a transfered into the governing 
differential equation. This is accomplished by the transformation 

g(x) = ah(t) with x = t 


(16) 

which leads to the initial value problem 



hh~ + t - L = 0 


(17) 

h(0) = 1, h'(0) = 0 

with 

| L (a , X) = a “2/3 


(18) 

(19) 

Let y(t) be the vector (y^t), y 2 (t)) t with y 2 = h(t) and 
Y2 = h'(t). Then, (17) and (18) can be written as the first order 

system 

dy/dt = F(t jy) 1 - 


(20) 

U 

O 

II 

/-N 

o 

v-/ 


(21) 

with F = (F 1> F 2 ) t and 



F l(t,y) = y 2 (t) 


(23) 

F 2(t,y) = (L - t)/y 1 (t). 


(24) 
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Thus, Fj and F 2 are continuous and bounded functions of t and y provided 
1 1 1 < °°, | y ^ | > 0, and ly^l < 00 . Further, with these restrictions on t and 
y, partial derivatives of F with respect to the components of y are 

continuous and bounded. Thus, F satisfies a Lipschitz condition on the set 
in (t.y)-space with |t| and |y 2 | bounded and |yj - l| < 1 - e with 

0 < e < 1. 

The Picard-Lindelof Theorem (Coddington and Levinson, 1955) now gives the 
existence of a unique continuous solution of (20) and (21), and hence (17) and 
(18) on the interval |t| < t Q with t 0 > 0. Further, for t in this 
interval, if h(t) is non-zero, then 

t 

h'(t) = / (L - s)/h(s)ds. (25) 

0 

In particular, if h(tQ— ) is non— zero, then h(tg) is non-zero and |h'(tg)| 
is bounded. The solution of (17) and (18) can thus be continuously extended 
past tQ , Two cases must now be considered* 

1. A continuous solution exists for all bounded t: Consider the range 
0 < t < L. Then, by (25), as h(0) is positive, both h(t) and h'(t) are 
positive, monotone increasing functions. Further, h'(t) is bounded above by 
L 2 /2 and has a local maximum at t = L. Beyond t = L, the factor L - s in 
(25) changes sign and becomes negative, and h'(t) becomes a monotone 

decreasing function which passes through zero at the point t* defined by 

L t* 

/ (L - s)/h(s)ds = / (s - L)/h(s)ds. 

0 L 


( 26 ) 
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For L<t<t*, h'(t) is positive and remains bounded above by L^/2. 
Consequently, h(t) is monotone increasing on this range and has a finite 
positive maximum value h at t . 

For t > t*, h'(t) is strictly negative and both h(t) and h'(t) are 

monotone decreasing. Two possibilities now exist. Either h(t) remains 

positive for all t > t*, or there is a finite T at which h(T) = 0. 

The first possibility leads to a contradiction. Suppose there is a 6 
such that 0 < 6 < h* < 00 and h(t) > 6 for all t > t*. Then, 

t t 

h'(t) = - / (s - L)/h(s)ds < -(1/6) J (s - L)ds. (27) 

* * 

t t 

This implies that as t tends to °°, h(t) remains positive but h'(t) tends 
to -°° , which is a contradiction. 

If the solution h(t) of the initial value problem can be continuously 
extended for all bounded t, there must thus be a finite T at which 

h(T) = 0. This, in turn, requires that if g(x) is the corresponding 

solution of the original initial value problem (10) and (15), then there is a 
finite X given by 

X = <x 2/3 T (28) 

at which g(X) = 0. As by (15) g'(0) = 0, g(x) will be a non-negative 

continuous solution of the boundary value problem (10) and (11) if X = 1 + X. 
Equation (28) now shows that one appropriate choice for the initial value 
a is 

a(X) = ((1 + X)/t) 3/2 . (29) 
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2. The solution h(t) cannot be continuously extended to all bounded 
values of t: In this case, there is a finite value t such that a 

continuous solution exists for 0 < t < t, but no continuous solution exists 
for t > x. 

To treat this case, we first note that t must be greater than t*. For 
if h(t) can be continued past tQ to some point t^ less than t*, then 
h(t) is monotone increasing, hCtj) is positive, and Ih'Ctj) | < l 2 /2. 

Consequently, h(t) can be continued past and in particular can be 

continued to t*. As h(t*) is positive and h'(t*) = 0, h(t) can now be 

continuously extended beyond t*. For t > t*, because t* > L, the term 
t - L in (17) is now strictly positive, so h'(t) will be strictly negative 
and both h(t) and h'(t) will be monotone decreasing. 

Two possibilities must now be considered. Either there is a T in the 
range t* < T < x such that h(T) = 0, or h(t) > 0 for t* < t < x. In the 
former case, we can proceed as before and define an appropriate a in terms 
of X and T by (29). Therefore, it is only necessary to consider the 

second alternative. 

Suppose now that h(t) is positive for t < x. Then, it is clear that 
the limit of h(t) as t tends to x exists and hence h( x) is well 

defined. Further, if h(x) > 0, then |h'(x)| is bounded. Under these 

circumstances, a continuous solution h(t) would exist in a neighborhood of 
x, and in particular, for some t > x. This would contradict the definition 
of x. Consequently, we must have h(x) = 0. An appropriate initial value 
o in terms of X and x is now defined by (29) with T replaced by x. 
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The theorem is now established. 

* 

The above proof shows that for 0 < X < X there is at least one 

positive initial value o which leads to a solution of the boundary value 

problem. In the case when X = 0, Callegari and Friedman showed that there 

was in fact exactly one positive a; i.e., the solution of their boundary 

value problem is unique. Their uniqueness proof is critically dependent, 

however, on the fact that the term 2uk(u) in their equation, which plays the 

role of the x - X term in equation (10), is non-negative over the full range 

of x values. By contrast, x - X in (10) clearly violates this property. 

* 

Indeed, as will be seen in the next section, for 0 < X < X , there are two 
positive values of X which lead to solutions of the boundary value 
problem. Hence, for this range of X, solutions of the boundary value problem 
(10) and (11) will exist but will not be unique. 


4* Non-uniqueness of Solutions 

To examine the dependence of solutions on X and a, the initial value 
problem 

f'~ + ff~ = 0 


f (0) = 0, f'(0) = -X, f~(0) = a 


(30) 


was solved numerically. This problem in the original similarity variables 
corresponds to the Crocco variable initial value problem (10) and (11). In 

this formulation, it is clear that the initial value a, physically, is 
related to the skin friction. 
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The initial value problem (30) was first solved for fixed X to determine 
the limiting value f"(°°) as a function of a. Results of the calculations 
are shown in Figure 1 . For values of a at which the curves cross the 
horizontal line f^(°°) = 1, a solution of the boundary value problem (1), (8), 
and (9) is obtained. 

Consider first the curve for X = 0, which is the usual Blasius problem 
considered by Calligari and Friedman. In this case, f"(°°) is a monotone 
increasing function of a which crosses the line f'(°°) = 1 at 
otg = .46960... . This value of c*q is well known from the homotopy 
technique. 

For X > 0, the curves in Figure 1 have a local minimum value of f'(°°) 

which will be denoted by M(X). Calculations show that M(X) is an 

^ & 

.increasing function of X. When X = X , M(X) = 1, and there will thus be 
exactly one value of a, a say, such that a solution of the initial value 


problem (30) is also a solution of the boundary value problem (1), (8), and 


(9). For X>X,M(X)>1 and there will thus be no solutions of the 

boundary value problem. When 0<X<X,M(X)<1, and hence there will be 

two values of a which give solutions of the boundary value problem, i.e., 

non-uniqueness. These values lie in the ranges 0 < < a and 

* 

a < a 2 < a Q . 

To accurately determine X and a and refine the points at which the 
curves in Figure 1 cross the line f'(<») = 1, the basic initial value solver 
for fixed X was augmented by an iteration procedure on a which terminated 
when |f'(°°) - l| was less than 10“^. Figure 2 shows the resulting o's as a 
function of X. Critical values were found to be 
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X* = 0.3541078... and a* = 0.2180238... . 


(31) 


The non-uniqueness of solutions of the boundary value problem (1), (8), and 
(9) for 0 < X < X is clearly evident in Figure 2 • 
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1.3r 



f!"(0)=d 

Figure 1. The limiting value f* (°°) as a function of the initial 
value f"(0) = a for A= = 0, 0.1, 0.2, X , and 0.4. 

At the critical value X = X * , the minimum value 
M(A*) =1. c 
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